On diffraction and absorption of X-rays in perfect crystals 
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This work reports the incompatibility of the Stokes relation for time reversibility of wave's propa- 
gation and the usage of complex atomic form-factor in diffraction theories. The implications behind 
a diffraction theory that obeys the Stokes relation is discussed. 
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I. INTRODUCTION 

Interaction properties of X-rays with atoms are de- 
scribable by complex atomic scattering factors in the 
widely used format 



f = fo + f + *f" 



(1) 



Atomic resonance theoriesi*2iMi& provide the physical 
bases for the /' and /" corrections of fo, the scatter- 
ing factor for incident photons with energies away from 
the absorption edges of the atom. 

In an aggregate of atoms where the scattering ampli- 
tudes of the individual atoms have random phases, the 
material's index of refraction is given by 



2wV 2-*> 



(fo + f + if"), 



(2) 



n runs over all atoms in the volume V, A is the X-ray 
wavelength, and r e = 2.818 x 1CP 5 A is the classical elec- 
tron radius. 

(1/^0 En(/o + /')n provides the effective electron den- 
sity, and it is responsible for the difference from unit of 
the real part of the index of refraction, and then, for the 
effect of refraction. The imaginary part of n r , propor- 
tional to J2 n f'ni i s responsible for the intensity attenua- 
tion of the travelling X-ray waves in the medium. 

In crystalline materials, the scattering and absorption 
properties of the atoms are summarized in the structure 
factor 



•F h = Y,(fo + f' + */")» exp(+27riff • r„) 

n 

= \F H \e ia », (3) 

which is obtained as the scattering amplitude by the crys- 
tal's unit cell for a given reflection H, whose diffraction 
vector is H . Atomic positions in the unit cell are speci- 
fied by r„. 

The above structure factor expression, Eq. is well- 
accepted since it is able to provide an explanation for 
different intensities of the H and H reflectionSfSi^ i.e. for 
the breaking of the Friedel's LawiS It is possible because 



in noncentrosymmetric crystals \!Fh\ ^ \J~ff \ owing to the 
/" scattering amplitudes, commonly called anomalous or 
resonant scattering amplitudes. 

There are several structure phasing methods currently 
in use based on the anomalous scattering phenomenon, 
such as multi-wavelength anomalous dispersion or iso- 
morphous replacements of anomalous scatters (see for 
reviews Refs. El El and These methods have al- 
lowed X-ray crystallographers to solve thousands of pro- 
tein structures, which is a strong evidence that complex 
atomic scattering factors properly describe the interac- 
tion physics of the X-ray waves with the atoms. More- 
over, near absorption edges in cases where fo + f = 0, 
the dynamical diffraction of only the resonant scattering 
term if" has been theoretically predicted and demon- 
strated experimentally«i2ii2iiii 

The initial motivation for this work was related to 
energy conservation in approximated solutions of multi- 
beam scattering problems^ when applying the Stokes 
relation of thin film optica^ to each lattice plane of a sin- 
gle crystal. It has led to an incompatibility with the us- 
age of complex atomic scattering factors, e.g. Eq. QJ, in 
diffraction theories that exists even in two-beam diffrac- 
tion cases, as further discussed here. 

An apparent incompatibility arrives because the reflec- 
tion and transmission coefficients derived from Stokes 
relation do not allow absorption to be treated as 
anomalous-resonant scattering amplitudes. But, on the 
other hand, these coefficients preserve the sum of proba- 
bilities for reflection, transmission and absorption of pho- 
tons at each lattice plane, and they provide an accurate 
description of the X-ray diffraction in perfect crystals 
with thickness varying since the kinematical regime of 
diffraction (few lattice planes) to the dynamical one (very 
thick crystals). Then, the point of conflict relies on how 
these coefficients that provide a good description of the 
diffraction phenomenon can not be in agreement with 
atomic resonant theories. 

To properly explain the implications behind a diffrac- 
tion theory that obeys the Stokes relation and to discuss 
its consequence on X-ray crystallography, this article has 
been divided into three Sections (II, III, and IV): 

Section II: Derivation of reflection and transmission 
coefficients from the Stokes relation for a general quan- 
tum particle. Recursive equations for calculating the re- 
flectivity of absorbing crystals with a finite number of 



lattice planes are also obtained. These recursive equa- 
tions are extremely simple and, even thus, capable to de- 
scribe important features of the diffraction phenomenon 
such as kinematical diffraction, primary extinction, in- 
trinsic width of Bragg reflections in either kinematical 
or dynamical diffraction regimes, and phase shift of the 
diffracted waves across the reflection domain. 

Section III: Application of the recursive equations, 
derived in Section II, to the particular case of X-ray 
diffraction. It emphasizes the formation of standing 
waves as necessary to explain the intensity differences of 
Friedel reflections in noncentrosymmetric crystals when 
the structure factor is calculated without the resonant 
term if" . 

For a given example, the anomalous signal (intensity 
difference from a pair of Friedel reflections) calculated 
here as an effect of absorption modulation by standing 
wave o 18 ' 19 i 20 ' 21 is compared to that obtained by the usual 
expression of the structure factor, as well as to experi- 
mental data. This comparison allows a clear distinction 
between absorption and interference effects owing to the 
if" term. It also demonstrates that the majority of the 
anomalous signal amplitude has its source on absorption 
modulation and therefore it should be a function of the 
thickness when the crystal is very small (below the pri- 
mary extinction length). Implication on the accuracy of 
phasing methods based on anomalous dispersion is dis- 
cussed. 

Section IV: A discussion is given on a possible physi- 
cal explanation for the apparent incompatibility of the 
Stokes relation [e.g. Eq. |0J] and the usage of the 
anomalous-resonant term in the expression of the struc- 
ture factor, i.e. the usage of if" in Eq. (J3J. The central 
point of discussion relays on causality, collapse of quan- 
tum probability amplitudes for photons, and stimulated 
photon emission by excited atoms. 



II. STOKES RELATION IN DIFFRACTION 
THEORIES 

The elastic scattering of quantum particles (photons, 
electrons, neutrons, ...) by a thin plane of matter can 
be described by introducing reflection R, and transmis- 
sion T, coefficients so that ipR = Ripo and ipT — Tipo 
are the probability amplitudes (or the vectorial electric 
field amplitudes in the case of photons) for the reflected 
and transmitted particles, respectively. When the prob- 
ability amplitude ipo of the incident particles is uniform 
over a significant area of the plane, regarding the dimen- 
sions of its building blocks, the particles are either spec- 
ular reflected or transmitted along the same direction of 
the incident ones, as for source Si in Fig. 1. |i?.| 2 and 
|T| 2 are therefore the reflection and transmission proba- 
bilities of detecting the particles along these directions, 
as indicated in Fig. 1 by the detectors Di or D2, re- 
spectively. For non-absorbing and non-dispersive planes 
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FIG. 1: Specular scattering of quantum particles (photon, 
electrons, neutrons, ...) from identical sources Si and S2, 
symmetrically displaced at both side of a thin plane of mat- 
ter. The reflected and transmitted probability amplitudes are 
given by ip R = Ripo, ipr = Tipo, tpn = Ripo, and ipf = Tip . 

\R\ 2 + \T\ 2 = 1, otherwise a = 1 - |i?| 2 - |T| 2 is the effec- 
tive absorption probability accounting for the probability 
of the particles to be in fact absorbed (photoelectron ab- 
sorption) or inelastic scattered at random directions. 

The reflection and transmission coefficients are in prin- 
ciple different, given by R and T, for incidence on the 
other side of the plane, as for source S2 in Fig. 1. Since 
the sources are identical and symmetrically displaced on 
each side of the plane, the sum of the counting rates 
from both detectors are equal to the total emission rate 
of the sources. It means that, for non-absorbing planes, 
\ipR + V'tI 2 + + iPt\ 2 = 2|-0o| 2 , and hence 

Rf* + R*f + RT* + R*T = 0. (4) 

It stipulates some phase relationships between the re- 
flected and transmitted amplitudes, as usually obtained 
for laser beam splitters 2 ^ or, equivalently, by the Stokes 
relation for time reversibility of wave's propagation^. 

Without loosing generality, the probability amplitude 
coefficients can be written as 

R = ±i\R\ e ii5+ *\ T= \T\e ilp , 

R = ±i\R\ e'G+ri, and T=\f\e^ (5) 

where according to Eq. (@J, 

|i?||T|sin<5+ |i?||T|sin<5 = 0. (6) 

ip and Cp are the phase delays across the plane, whose 
thickness d is comparable to the particle's wavelength A. 
Since the incidence angle 9 is the same on both sides 

(p = (p = — — dsinO . (7) 



<5±90° and <5±90° are the amount by which the phases 
of the reflected particles differ from the phases of the 
transmitted ones. If S = 6 = 0, the phase of the reflected 
particles are 90° shifted with respect to the transmitted 



ones 



1,23.24 



A. Reflection and transmission coefficients for 
N-plane crystals 

In crystals, each individual scattering plane stands for 
an element of periodicity composed, in general, of several 
atomic layers that repeat themselves with period d in 
order to build up the crystal. The scattering properties of 
each element of periodicity (lattice plane) , when specified 
by reflection and transmission coefficients, Eq. (J3J), allow 
to calculated the Rn, Tjv, Rn, and Tjy coefficients for 
crystals with N planes in reflection geometry. 

To carry out this calculation, let first analyze the two- 
plane case. In such case, the particles can suffer several 
reflections in-between the planes before exiting, as de- 
picted in Fig. 2. However, even when there is no empty 
space h between the planes (h — ► 0), the contributions of 
the several consecutive reflections have to be taken into 
account. It provides that for N = 2 




T\ t\ t\ t\ 

FIG. 2: Reflection and transmission coefficients of two planes. 
Several consecutive reflections occur in-between the planes, as 
in the resonant cavity of a Fabry-Perot interferometer. They 
have to be taken into account even when there is no empty 
space h, between the planes, i.e. when h — > 0. R, T, R, and 
T are the amplitude coefficients of each individual plane, as 
given in Eq. @. 



are the amplitude coefficients for A-planes in a crystal 
of thickness Nd. 



B. Absorption 

The absorption probability 



oo 

i? 2 = R + TRT + TRRRf + ... = R + TRT ^ (RR) n 

n=0 

(8a) 

and 

oo 

T 2 = TT + TRRT + TRRRRT + ... = TT^(RR) n 

71=0 

(8b) 

for incidence from the top-left, as in Fig. 2. In the 
other symmetrical situation, which corresponds to inci- 
dence from the bottom-left, e.g. source S2 in Fig. 1, the 
i?2 and X2 coefficients are obtained by analogous pro- 
cedure. Note that these coefficients are very similar to 
those in the Airy's formula of the Fabry-Perot interfer- 
ometer, and that their expressions can be simplified by 

usin g Y,n=o 7 - n = !/(! - z ) for \A < L 

Although it would be possible to add plane-by-plane 

with the above procedure, a fast recursive formula to go 
from very thin to semi-infinity crystals (of infinity thick- 
ness) is obtained by building the crystals in geometrical 
progression of N = 2™ planes (n — 1,2,... ). Since the 
equations to calculate the coefficients of N planes from 
those coefficients of N/2 planes are identical to Eqs. ijHJ), 



a = fid/ sin9 (10) 

is given in terms of the linear absorption coefficient fi of 
the material. It does reduce the transmission probability 
per lattice plane according to 

\T\ 2 = (1 - |i?| 2 ) e- a ~ (1 - |i?| 2 )(l - a) ~ 1 - |i?| 2 - a 

(11) 

since a\R\ 2 < |i?| 2 + a < 1. 

C. Intrinsic width and phase shift in rocking curves 

As a function of the incidence angle 8 — the variable 
angle in rocking curves — the reflection coefficient of N 
planes, 

Rn(6) = \R N (6)\e^W (12) 

determines the intrinsic profile of the diffraction peak, 
i.e. the intensity reflectivity curve \Rn{0)\ 2 , as well as 
the phase 



Rn 
Rn 



and 



= 1 



1 - Rn/2Rn/2 



1 — RNI2R 



N/2^N/2 



N/2 
N/2 



T 2 
± N/2 

± N/2 



(9a) 



(9b) 



*(0) = 12(0) + 5, 



(13) 



of the reflected particles with respect to the incident ones. 
0,(9) stands for the phase dependence with the incidence 
angle, and 5 is determined by the internal structure of 
the lattice planes. 



The center of the reflectivity curve is observed at the 
Bragg angle 8b- At this angle, the phase of TT in 
Eq. H8a|l is an integer number of 27r, i.e. 2ip = —2mn 
(to = ±1, ±2, ...) and hence, for the first order reflection 
(to = 1) we have that sin 9b = X/2d. 

A few examples on the behavior of the intensity re- 
flectivity curve \Rn(0)\ 2 as a function of crystal thick- 
ness, reflection and absorption probabilities are shown 
in Fig. 3. The phase shift fl(8) across the total reflec- 
tion domain of the rocking curve is also shown (inset of 
Fig. 3). In all calculated curves, \R\ = \R\ and \T\ = \ f\. 

Absorption is observed to symmetrically round (re- 
garding the center of the peak) the flat-top of the total 
reflection domain, and to reduce both the integrated in- 
tensity (area of the rocking curve) as well as its maximum 
value. 

The analysis of the reflectivity profiles for thick crystals 
leads to an empirical equation for the full width of the 
half maximum (FWHM), 



W=- x /\R\\R\tan0 B . 



(14) 



For a given Bragg angle, the width of the rocking curve is 
essentially determined by the particles specular reflection 
probability i?| 2 , at each lattice plane. When the interac- 
tion physics of the particles with the planes are known, 
the FWHM can be calculated by means of more specific 
parameters such as electronic density, atomic numbers 
(neutron scattering), X-ray polarization, particle energy, 
etc. On the other hand, the reflection probability can 
be estimated only from the values of W and Ob ■ For in- 
stance, the silicon 111 reflection at about 8KeV presents 
W ~ 10" and an incidence angle close to 14°, providing 
\R\ 2 ~ 8.5 x 1(T 8 . 



III. X-RAY DIFFRACTION 

In the case of X-ray diffraction, the R and R reflection 
coefficients of a single lattice plane should not bear the 
resonant term if" . Otherwise, in general, the relation- 
ship in Eq. © will not be fulfilled in crystals without 
a center of symmetry (see Section IV for more details). 
To accomplish the condition imposed by this relation- 
ship, the scattering and absorption of X-rays by individ- 
ual atoms are treated here as follow. 

The atomic scattering amplitude taking place in inter- 
ference phenomena is given by 



a s (2d)=r e X(fo + f')\C\, 
while the atomic absorption cross section 
a A = 2r e Xf" 



(15) 



(16) 



accounts for photoabsorption. \C\ = 1 or cos 29 stands 
for the a or tt polarization components, respectively, /o 
is the main scattering factor through an angle 29 from 
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FIG. 3: Intensity reflectivity curves \Rn(0)\ 2 in crystals of 
thickness Nd, according to Eqs. ©. A9 = 6 — 8b, A = 
1.54A, and d = 3.14A in all of the following cases: (1) [N, 
\R\ 2 , nicm,- 1 )] = [2048, 16 x 10~ 8 , 0.01] (pink curve); (2) 
[4096, 16 x 10~ 8 , 0.01] (blue curve); (3) [00, 16 x 10~ 8 , 0.01] 
(red curve); (4) [00, 64 x 10~ 8 , 0.01] (green curve); (5) [00, 
64 x 10~ 8 , 200] (orange curve); (6) [00, 64 x 10~ 8 , 8000] (black 
curve); and (7) [4096, 64x 10~ 8 , 400], for which both \Rn{6)\ 2 
(dashed curve) and \Tn(6)\ 2 (black-line curve) are shown at 
the top-left inset. For cases (3), (4), (5) and (6), the phases 
of Rn(8) are given in the top-right inset; they correspond to 
Q(8) since 8 — in Eq. I13H was used in all simulations. 



the incident beam direction. More details on /o, /', and 
/" values are available on the International Tables for 
Crystallography^ 

A single layer of atoms, all of the same kind, with M 
atoms per unit area scatters an incident plane wave of 
amplitude Eq according tc 23,24 



Esl — 



. M 

1 

sin 9 



<js{29)E = RslEq. 



(17) 



\Rsl\ 2 is the probability to measure the specular scat- 
tering of photons by this single layer of atoms, while the 
absorption probability is given by 



M 



asL 



sin t 



-,0-A- 



(18) 



To generalize the reflection coefficient Rsl, from the 
single layer to elements of periodicity composed of sev- 
eral atomic layers, as illustrated in Fig. 4, the number of 
atoms per unit area is written in terms of the crystal's 
unit cell, so that 



Rsl = -i- 



M 
sin 9 



o-s 



R 



V c sin 9 



-2-KiH- 



. ( 19) 

n runs over all atoms in the unit cell, whose positions 
and scattering amplitudes are given by r n and <Ts(n), 
respectively. H is the diffraction vector of the chosen 



reflection H. V c is the unit cell volume and d = 1/\H\ 
as usual. The exponents, or geometrical phase factors, 
are necessary to account for the phases of the photons 
coherently scattered at the atomic sites, and its " - " 
(minus) signal will be justified later on. 

By replacing Eq. (|15J) into Eq. (|19|l and defining 



F h = + /')nexp(-27rifr • r n ) = \F H \e l5H (20) 

71 

as the structure factor of reflection H, the reflection co- 
efficients for individual lattice planes arc 



j j jj * * * * 
'> jjj j j j j j j j j 

> ** * j * j j 

mi ******** mm 

r / j * jj * * * j 

/ * **> ** j a ^ 



**» *mm ** *j j jj 

jj jj j j j j 
******* j jj * 
* ** j j j j .j 



H 



r e X\C\d 
V c sin(9 



(21) 



where -Fg 



stands for the structure factor of reflection H 
with diffraction vector — H . Temperature factors are im- 
plicit in the / and /' values. The photoabsorption prob- 
ability by the lattice plane is obtained with the same 
procedure, i.e. 



FIG. 4: Each crystal reflection stands for an element of peri- 
odicity composed in general of several atomic layers. Regard- 
less the chosen origin, the r n atomic positions in the elements 
of periodicity (lattice planes) are specified from bottom to 
top, in the same sense of the diffraction vector H . This is the 
usual convention in X-ray crystallography. 



A. Absorption modulation by standing waves 



asL 



M 
sin 9 



OA 



V c sin(9 



]cr A (n). 



(22) 



IS 



No geometrical phase factor is necessary since oa{ji 
not a scattering amplitude. 

To make certain that the above results are correct, 
they can be compared with those already known from X- 
ray diffraction theories. The width of the total reflection 
domain, characterized by top- hat shape of the reflectivity 
curves in Fig. 3, also known as the Darwin width^2SiSZiS£ 
is obtained by replacing Eq. l|2T|l into Eq. (|T4"f> . and then, 
the FWHM is given in terms of the same parameters used 
to calculate the Darwin width Wd, he. 



ff , (23, 

3 V c cos 3 

see Ref . |2t| for a comparable expression of Wd ■ 

The phase shift of 180° across the rocking curve, Q(0) 
in the inset of Fig. 3, agrees with that calculated by the 
dynamical diffraction theory, which is known as the dy- 
namical phase shift 30,31 . Only the exact center of the 
reflectivity curves, at 9b, are slightly different since A 
has not been corrected for the effect of refraction. 

Eqs. I |10| l. (|16|l and 1221) can be combined to provide 
the linear absorption coefficient as 



In order to fulfill the condition imposed by the Stoke 
relation, e.g. Eq. (|34|l . structure factors were calcu- 
lated without the resonant terms if^. Consequently, 
another explanation for the observed anomalous signals 
from Friedel pairs is needed; one that is not based on the 
differences of the square modulii of the structure factors. 

One possible explanation comes from that old interpre- 
tation on the contrast of Kossel lines provided by LaueiS 
and on Borrmann anomalous transmissiorii2i22i2i were 
the photoabsorption is proportional to, or modulated by, 
the intensity of the standing wave field formed in crystals 
undergoing diffraction. 

However, to a general explanation of the anomalous 
signal based on standing waves, there is an important 
difference regarding the current understanding on this 
subject jSL^iiSSiSM the positioning of the standing wave 
field must be independent of the absorbency of the 
crystals 3 ^, as well as invariant regarding arbitrary choice 
of origin for the atomic positions^ A detail discussion on 
the dependence of the standing wave field with the phase 
of the structure factor has been carried out by Bedzyk & 
Mater lik 34 . But, on its usual expression the phase an of 
the structure factor in Eq. Q depends on absorption. 

To demonstrate that absorption modulation by stand- 
ing waves (AMSW) can provide a significant contribution 
to the anomalous signal, the incident Ei nc and reflected 
Efief waves as a function of depth h are approximated 
to 



(24) 



E 



Inc. 



= E (h)v e l{u)t - k °- r) 



(25) 



and 



which is the usual expression for calculating [i of a crystal 
in terms of its unit cell, see for instance Ref. 0. 



E Ref = R N (9)E (h)v e*M-<-) (26 ) 
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FIG. 5: Standing wave positioning must be invariant regard- 
ing an arbitrary choice of origin—. Such requirement is ac- 
complished in the present approach when the structure factor 
is calculated with the " - " signal in the phase factor. For a lat- 
tice plane made of a single atomic layer A, the structure factor 
phase is given by (a) Sh = —2irzo/d or (b) Sh = +27r«o/d. 
As the crystal is rocked through reflection H, the standing 
waves are formed at the positions given by dark lines, move 
across the gray shaded areas, and vanish at the dashed-line 
positions. In (a) the scanned portions of the lattice planes 
(gray areas) are invariants regardless the value of Zo, but the 
same invariance is not observed in (b). 



where fco and k are the wavevectors of the waves, while i>o 
and v are their oscillation directions so that, |C| = i) ■ i>o 
for both a and ir components of polarization. A com- 
parable reflection coefficient at all depth is assumed, i.e. 
Rn(9) do not depend on h. All dependence with depth is 
accounted by EqQi — > oo) — > 0, which is a smooth func- 
tion with very small variation over the lattice period. 

The time average intensity of the standing waves is 
then 



E 



Ref. 



\E (h)\ 2 {l + \R N (6)\ 2 



(27) 

■2|i?jv(0)||C|cos$} 



where $ = 2nH ■ r + *(0), 2itH = k - k , and 
2itH ■ r = 2-nz / d. The node positioning along the lattice 
plane normal direction z, which is parallel to the diffrac- 
tion vector H = [l/d)z, is obtained from the condition 
$ = 27T, so that 



. . (d(l/2-<W27r) forn = 7T, 
Z ^ = ^-^^={ d \l-sjl) for O = 0. 

(28) 

Here, the positioning of the standing waves are already 
independent of absorption since Sh came from the defi- 
nition in Eq. H2()Jl : and the condition of invariance with 
an arbitrary choice of origin is fulfilled by the " - " signal 
of the phase factor in Eq. (|19fl . It can be demonstrated 
by the following example. 

For a given reflection, whose element of periodicity is 
made of a single atomic layer A, as shown in Fig. 5, the 
origin is taken arbitrarily so that it does not necessarily 
coincident with the atomic layer, i.e. zq =/= 0. When the 
structure factor phase is Sh = — 2Trzo/d 1 according to 
Eq. (|28fl. the node position moves from z no d e = d/2 + z 
(ft = 7r) to z noc ie = d + z (SI = 0) across the rocking 




H 



B A 



FIG. 6: In noncentrosymmetric crystals, A/B, the standing 
waves glide over regions (gray areas) with different electron 
densities during the rocking-curves of the (a) H and (b) H 
reflections. In centrosymmetric crystals, A = B, identical 
regions are scanned by the nodes 



curve, as represented by the gray shaded area in Fig. 5(a). 
It is an invariant behavior since the nodes will always 
scan the same portion (gray area) of the lattice planes 
regardless the value of zq. On the other hand, the same 
invariant behavior does not occur if the geometrical phase 
factor is calculated with the " + " (plus) signal, i.e. Sh = 
+2irzo/d, as shown in Fig. 5(b). 

One may note that the gliding sense of the nodes, in 
the sense of the diffracting vector in Fig. (5), is opposite 
to that reported on most articles on this matter. The 
gliding sense is determined by the rotation sense of the 
dynamical phase shift 0,(6). A clockwise rotation sense 
of SI (from 7r to 0) is that predicted by the dynamical 
theory (e.g. Ref. l30f) . and it is the same of that obtained 
here for the Rn(&) coefficient (see inset of Fig. 3). Then, 
the gliding sense of the nodes here should be the same of 
that reported by other authors. 

Assumption of an opposite gliding sense, in the sense of 
H , would be necessary to math fluorescence results mon- 
itored during a rocking curve with the integrated inten- 
sities of a Friedel pair calculated according to Th (with 
the " + " signal in the phase factor) £L 

For the invariance of the wavefield regarding the choice 
of origin, either a clockwise or counterclockwise rotation 
sense of SI would require the " - " signal in the phase fac- 
tor. And, to avoid any confusion regarding the gliding 
sense of the nodes in the remaining parts of this article, 
the geometrical phase factor of Th is henceforth calcu- 
lated with the " - " signal . 

Since the nodes glide half a lattice plane distance as 
the crystal is rocked across the reflection domain, the 
portion of the lattice plane scanned by the nodes, e.g. 
the gray areas in Fig. 5, can be different for the H and H 
reflections depending on the symmetry of the crystalline 
structure. 

In centrosymmetric crystals the structure factor phase 
when calculated with the " - " signal in the phase 
factor — set the nodes to scan the lattice planes from the 
symmetry center with lower electron density to the one 
with higher electron density^ It implies that, equivalent 
portions of the lattice planes are scanned during rocking- 
curves of Friedel related reflections and, consequently, 
the effective photoabsorptions are exactly the same for 
both reflections. It is illustratively shown in Fig. 6 when 



TABLE I: Intensity variation of hkl /hkl Friedel reflections 
in a GaSb crystal at photon energy of 10.0 KeV. Qf and Qa 
stand for theoretical anomalous signals when the integrated 
reflectivities are taken as proportional either to |-7"if| 2 [Eq. 
with the " - " signal in the phase factor] or to Ah [Eq. I|33[l ] . 
respectively, r = Qf/Qa and r* regards this ration at 10.4 
KeV (above the Ga edge at 10.37KeV). Crystal dimension 
larger than 20/^m in all directions (see Fig. 8). 



hkl 


Qf (%) 




Qa (%) 


Ah /Ah 


r 


r 


111 


-5.817 


0.890 
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1.48 


113 


+6.845 
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1.102 


1.41 


1.52 


133 
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-5.255 


0.900 


1.42 


1.55 


115 


-7.906 


0.853 


-5.519 


0.895 


1.43 


1.58 


135 


+8.202 


1.179 


+5.688 


1.121 


1.44 


1.61 


117 


+8.450 


1.184 


+5.786 


1.123 


1.46 


1.67 


555 


-8.209 


0.848 


-5.531 


0.895 


1.48 


1.76 



assuming identical A and B layers of atoms. 

On the other hand, in noncentrosymmetric crystals, 
the scanned portions of the lattice planes for the H and 
H reflections may not be exactly the same and then, the 
effective photoabsorptions for some Friedel pairs can be 
different. This situation can also be illustrated by the 
example in Fig. 6 assuming now that the A and B layers 
of atoms have different absorption properties. 

To calculate the effective absorption coefficient 
u, e ff{0), as a function of the rocking angle, the normal- 
ized intensity of the standing waves 



Isw 



lIo I sw(z)dz 
2\R N (d)\\C\ 
+ l + \R N (9)\i 



cos(27ri? • r + *)(29) 



is used as a weight function 2 — for the atomic absorption 
cross sections in Eq. (I24[l . which provides 



Me//(f) 



2r e A 

Vc. 



(30) 



By implementing the absorption probability 
a = fj, e ff(0)d/ sin 9b, the reflectivity curves show 
their well-known asymmetric profiles with lower inten- 
sity at the right-hand shoulder 

For the sake of comparison, let define the amplitude of 
the anomalous signal from reflections H and H by 



Qf = 



\Fh\ 



\^h\ 2 



as in structure determination methods, and 



Qa = 



A H -Aa 



H 



A 



H 



An 



(31) 



(32) 
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FIG. 7: Anomalous signals from (a) 117/117 and (b) 115/115 
Friedel pairs in a GaSb crystal as a function of photon en- 
ergy across the Ga edge at 10.37 KeV. Experimental values 
were determined by measuring these asymmetric reflections 
at either low (squares) and high (circles) incidence angles. 
Theoretical Q f (solid lines) and Qa (dashed lines) curves are 
provided by Eqs. I3H and 1321 . respectively. Sample: com- 
mercial GaSb(OOl) wafer, 500^m thick, and polished both 
sides. X-ray source: synchrotron radiation, a polarization, 
and double-bounce Si(lll) monochromator. 



where 



A 



H 



\R N (0)\ 2 dO. 



(33) 



H 



Rn(6) is given in Eq. i|12|) for both reflections but cal- 
culated for different absorption probabilities as for the 
theoretical values in Table I and in Fig. 7. 

Since the anomalous signal in the AMSW hypothesis 
is related to the formation of standing waves, Qa is a 
function of the crystal thickness. Fig. 8 shows both the 
integrated reflectivity and Qa as a function of thickness 
for a few Friedel reflections in the GaSb crystal. It varies 
since a very small crystal undergoing kinematical diffrac- 
tion where the intensity increases with thickness, until 
the dynamical diffraction regime dominated by primary 
extinction where, for an absorbing crystal, the intensity 
no longer depends on thickness. 

Essentially, the observed disagreement between the Qa 
and Qp values in Table I, or in Fig. 7, arrives because in 
the former treatment interference/diffraction effects in- 
volving the if" resonant amplitudes are not taken into 
account. Below the Ga edge only fg b = 4.09 has a sig- 
nificant value regarding f^ a = 0.52, while above this 
edge they have nearly the same value, f'g b — 3.833 and 
f f Q = 3.87, and then the disagreement is more evident 
and beyond experimental errors. It demonstrates that 
the resonant amplitudes are taking part of the diffrac- 
tion phenomenon. 

On the other hand, the comparison in Fig. 7 also in- 
dicates that more than 50% of the anomalous signal am- 



8 




Thickness (um) 

FIG. 8: (a) Mean integrated reflectivity / = C(A H +Ag)/2, 
and (b) anomalous-signal amplitude Q — \Qa\, as a func- 
tion of crystal thickness. Friedel pairs: 111 /111 (1), 113/113 
(2), 133/133 (3), and 155/155 (4). Normalization constant 
C _1 = 85/i radians. Crystal: GaSb. Photon energy: 10 
KeV. 



plitude Q relies exclusively on the effects of absorption 
modulation, which implies that Q may depend on the 
crystal dimensions. 

As pointed out in Fig. 8, the relative behavior of Q as a 
function of thickness regarding strong {Q strong) an d weak 
(Qweak) reflections can lead to a procedure to check the 
AMSW effect on experimental Q values, i.e. kinematical 
diffraction: Q strong > Qweak', dynamical diffraction: 

Qstrong < Q weak ■ 

In low-absorbing crystals the primary extinction length 
can be of the order of hundreds of microns. Correcting 
the experimental anomalous-signal by a thickness related 
factor might be necessary to improved the accuracy of 
phasing methods in such small crystals. 



IV. ATOMIC RESONANCES AND 
INTERFERENCE EFFECTS 

The incompatibility of the Stokes relation, Eq. 
with the usage of if" in the expression of the structure 
factor can be demonstrated by replacing Eqs. JSj and 
(I21|l into Eq. iJBJ, which provides 

\J~h \ sin an + sinag = (34) 

since |T| and \T\ differ from unit by less than 10~ 6 . 

In noncentrosymmetric crystals, this equality is not ac- 
complished for Friedel reflections where ^ \J~h\i as 
can be easily checked by just replacing a few values. It 
has been very difficult to understand the physical reason 
behind this incompatibility. Does it prove that the res- 
onant amplitude if" should not appear in the calculus 
of the structure factors? But then, how to explain all 
the evidences mentioned at the Introduction section in 
favor of using the complex atomic scattering factor / in 




FIG. 9: (a) Interference of probability amplitudes of all pos- 
sible photon trajectories from the source (S) to the detector 
(D) providing the specular reflection (dark line trajectory) by 
an uniform plane of atoms, see for instance Ref. I37L chapter 
2. (b) When an atom of the plane is excited by the incident 
photon, all other possible trajectories collapse into the sin- 
gle trajectory that goes from S to the excited atom. Then, 
the pattern of interference providing the specular reflection is 
destroyed, and any further atomic scattering, generated as a 
consequence of the absorbed photon energy, should not inter- 
fere with the scattering from other atoms. 



diffraction theories, as well as the experimental results in 
Fig. 7? 

In the calculation of anomalous X-ray scattering fac- 
tors, "the scattering amplitude for light by a bound elec- 
tron is proportional to /. The value of / is computed as 
a sum over all intermediate electron states except those 
states occupied by other atomic electrons, and a further 
sum must be taken over all the electrons of the atom. 
When the set of intermediate states is complete, the sums 
provide /o the main part of the scattering factor for X- 
rays." 

The above paragraph summarize the initial assump- 
tions made by Cromer & Liberman* in their paper on 
anomalous X-ray scattering. It clearly demonstrates that 
the anomalous-resonant scattering with amplitude if" is 
a privilege of excited atoms. 

Since the energy exchange between atoms and the elec- 
tromagnetic fields is quantified, excited atoms inevitably 
imply in previous photoabsorption, which requires the 
photons to manifest themselves as particles. Conse- 
quently, the probability to find the absorbed photon in 
the entire space occupied the electric field has instanta- 
neously collapsed at the position of the excited atom, as 
pictorially shown in Fig. 8(b). Since the quantum proba- 
bility amplitude has been collapsed at one atom, the reso- 
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nant amplitudes from this atom should not interfere with 
the scattering amplitudes from other atoms, unless some 
stimulated emission mechanisms are involved. Therefore, 
would not be plausible to assume that the anomalous- 
resonant scattering with amplitude if" , when occurring 
only as a consequence of photoabsorption at random in- 
stants of time, should not take part in interference phe- 
nomena such as X-ray diffraction? 

This is one possible explanation capable to conciliate 
both the Stoke relation in Eq. (|34[) and atomic reso- 
nant theories predicting complex atomic scattering fac- 
tors, such as / = fo + f'+if"- To be in agreement to this 
explanation, the if" term has to be omitted in the struc- 
ture factor expression, as for instance in the calculation 
of the specular reflection coefficients given in Eq. I|21(l 
where only the atomic scattering amplitudes fo + f take 
part of interference effects^ Accounting for stimulated 
emissions in phase with the standing wave field may be 
attempted if the anomalous signal dependence with the 
crystal dimension is experimentally confirmed. 



curve 



13,23,26,27,28 



is obtained by a simple set of recursive 
equations valid for absorbing crystals of arbitrary thick- 
ness. However, this description is not conciliable with 
the usage of the anomalous-resonant term in the struc- 
ture factor expression. An analogous explanation of that 
given by Lau o 1 ?' 1 ? on the contrast of the Kossel lines is 
necessary to explain most of the intensity difference of 
Friedel reflections in noncentrosymmetric crystals. 
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